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In a hybrid Superconductor - Insulator - Normal metal tunnel junction biased just below the gap, 
the extraction of hot electrons out of the normal metal results in electronic cooling effect. The 
quasiparticles injected in the superconductor accumulate near the tunnel interface, thus increasing 
the effective superconductor temperature. We propose a simple model for the diffusion of excess 
quasiparticles in a superconducting strip with an additional trap junction. This diffusion model has 
a complete analytic solution, which depends on experimentally accessible parameters. We find that 
the accumulated quasiparticles near the junction reduce the efficiency of the device. This study 
is also relevant to more general situations making use of superconducting tunnel junctions, as low 
temperature detectors. 



In a Superconductor (S) - Insulator (I) - Normal metal 
(N) junction biased below the superconducting energy 
gap A, the quasiparticle tunnel current transfers heat 
from the normal metal to the superconductor. This re- 
sults in a significant cooling of the normal metal elec- 
tronic population^ In most cases, the cooling is less effi- 
cient than expected, in particular near the optimal bias 
just below the gap. It is often argued that a parasitic 
power could explain this behaviour. A more intrinsic 
explanation is local overheating of the superconducting 
electrode, as excess quasiparticles with an energy close 
to the gap accumulate near the tunnel contact and can 
tunnel back into the normal metal. The use of quasi- 
particle traps on the superconducting electrode has been 
found to improve the cooling performance^ A quantita- 
tive description of the overheating effects calls for a pre- 
cise analysis of the decay mechanisms of the excess quasi- 
particle population in the superconducting electrode. In 
this framework, a detailed theory of non-equilibrium phe- 
nomena in a micro-cooler with normal metal traps on the 
superconducting electrodes was recently developed:- 

This paper is an attempt to describe the low tem- 
perature diffusion of non-equilibrium quasiparticles in a 
superconducting strip coupled to a normal metal trap, 
in the presence of recombination and pair-breaking pro- 
cesses. Our model is an extension of the Rothwarf- Taylor 
equations^ for the relaxation of non-equilibrium quasi- 
particles. A normal metal trap junction provides an addi- 
tional channel for near-gap quasiparticle evacuation. Un- 
der realistic assumptions, the non-linear diffusion equa- 
tion can be solved analytically. The spatial profile of the 
quasiparticle population in the S strip is obtained. Near 
the tunnel junction, the quasiparticle density exceeds the 
thermal equilibrium density, even at low injection cur- 
rent. The heat transferred back to the normal electrode 
is found to reduce the effective cooling power of the tun- 
nel junction. We conclude by comparing our predictions 
to experimental results. 

The system we consider is the superconducting section 
of a low temperature N-I-S cooler (Fig. [lj. It is mod- 
eled by a semi-infinite line in contact with a dielectric 
substrate. The superconductor is covered by a second 
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FIG. 1: Schematic cut view of the system considered. The 
excess quasiparticles injected from the N-metal in the S elec- 
trode diffuse out of the junction region and relax by two pro- 
cesses: recombination and escape in the normal metal trap. 
The trap layer is supposed to overlap completely the super- 
conducting strip and both have infinite length. Both the sub- 
strate and the trap layer are at the bath temperature is To. 

normal metal N tra p separated by a tunnel barrier. A 
quasiparticle current is injected through a N-I-S junction 
at the origin of the x-axis. We assume that the layers 
thicknesses are small, so that we consider diffusion only 
along the x-axis. Practical micro-coolers are usually of a 
S-I-N-I-S geometry, i.e. they are symmetrically made of 
a double N-I-S junction. The geometry under study here 
thus corresponds to one half of a micro-cooling device. 

At zero current, the system is at equilibrium at the 
bath temperature To, assumed to be much smaller than 
the superconductor critical temperature. Upon increas- 
ing the current in the N-I-S cooler, the electron tempera- 
ture Tjv of the N-metal decreases while the quasiparticle 
population increases on the superconductor side. If the 
cooler junction is biased below the gap value, only the 
high energy tail of the normal metal electron popula- 
tion can tunnel. Accordingly, the distribution of injected 
quasiparticles is sharply peaked at the gap energy A. Its 
width of the order fcsTjv <C A (fee is the Boltzman con- 
stant) will be ignored in the following. 

The process of recombination in a non-equilibrium su- 
perconductor was discussed by Rothwarth and Taylor~ 
for thin film geometries. The decay of the quasiparti- 
cle density involves quasiparticle recombination retarded 
by phonon retrapping. In this process, two quasiparti- 
cles initially recombine to form a Cooper pair, resulting 
in the creation of a phonon of energy 2 A. This phonon 
can be subsequently reabsorbed by a Cooper pair, result- 
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ing in two new quasiparticles. The 2A energy leaves the 
superconductor when either the phonons or the quasipar- 
ticles escape to the bath. In our model, we introduce two 
additional contributions to the Rothwarth- Taylor equa- 
tions: the diffusion of quasiparticles in the superconduct- 
ing electrode and their absorption by the trap junction. 
In the superconducting strip, the system of equations for 
the density of quasiparticles per unit volume N qp and the 
density of 2A phonons A2A then reads: 
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Here D qp is the quasiparticle diffusion constant. It differs 
from the free electron diffusion constant because of the 
reduced group velocity of quasiparticles as compared to 
the Fermi velocity— There is no source term in Eq. 1 
and 2 as the source will be introduced as a boundary 
condition for the density gradient at coordinate x = 0. 

In Eq. 1, the first term in the right hand side is the 
recombination term. Its quadratic dependence indicates 
that the recombination is a two-particle process. The 
prefactor R is a constant depending on the electron- 
phonon interaction. The related recombination time 
T~, 1 = 2RN qp Q is strongly dependent on the temperature. 
At low temperature, it diverges as exp (A/To)£ The sec- 
ond term describes the quasiparticle creation by pair- 
breaking processes induced by a 2A phonon (two quasi- 
particles per phonon) . The related phonon pair-breaking 
time tb is weakly dependent on temperature, tr and tb 
are intrinsic characteristic times of the superconducting 
material ruled by the electron-phonon interaction^ 

Eq. 2 has similar terms for the reverse processes but 
without gradient term because the 2A phonon absorp- 
tion time tb is very short The low energy (compared to 
2 A) phonons are ignored since their interaction with the 
superconducting condensate is negligible in a BCS super- 
conductor at low temperature^ The trap film is assumed 
to be at the bath temperature. The escape terms on the 
right hand side of each equation describe respectively the 
quasiparticle escape by tunneling to the trap film at a 
rate t^ 1 and the phonon escape to the substrate at a 
rate r" 1 . The characteristic time To — R tra pe 2 N(EF)ds 
is proportional to the thickness ds of the S-film, the elec- 
tron density of states at the Fermi level N(Ep) and the 
specific resistance of the trap tunnel junction Rtrap- The 
escape time r 7 for 2A phonons to the substrate depends 
on the acoustic matching at the interface to the substrate 
and on the superconducting film thickness. 

With no injection current, Eq. 2 reduces to a simple re- 
lation between the equilibrium populations (labelled with 
a "0" index): RN qp0 = 2N 2A0 /t b . Following the BCS 
theory, N qp o decays exponentially at low temperature Tq 
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In the following, we will use the dimcnsionlcss relative 
density z(x) defined by: 



N qp {x) = N qp0 [l + z(x)} 



(4) 



Under steady state injection, Eq. 2 can be rewritten to 
express the density A^a as a function of z(x). Eq. 1 then 
gives the following second order differential equation for 
z(x): 
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where r e // = Tr(1 + t 7 /tb) is an effective quasipar- 
ticle recombination time^ In particular, the situation 
t 7 > tb corresponds to 2A phonons being re-absorbed 
before leaving the superconductor, which enhances the 
quasiparticle density. 

The boundary conditions at the two wire ends are: 
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The first one states that the excess quasiparticle density 
must decay to zero at the right end of S. The second one is 
the continuity condition for the quasiparticle flux at the 
injection point x = 0: the gradient of the quasiparticle 
density is equal to the rate of quasiparticle injection per 
unit area. Here A is the junction area and i is the bias 
current. The injected electrons create quasiparticle exci- 
tations with energy about A per electron, and no charge, 
the electrical continuity being insured by the supercur- 
rent in S. Charge imbalance effects can be ignored here 
since the N-I-S junction is biased near the gap edge.— 

The solution of Eq. [5] gives the analytic expression of 
z as a function of the x-coordinate: 
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It includes an integration constant a which depends on 
the injection current i. It is obtained by substituting 
z(x = 0) in the following equation: 
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Here we have introduced the parameter a = 1 +T e ///To, 
which is the enhancement ratio of the quasiparticle de- 
cay rate due to the presence of the trap junction. The 
quasiparticle diffusion length A = yl D qp T e f / ja combines 
the two mechanisms of absorption/recombination with 
phonons and trapping in the normal metal trap. When 
trapping is dominant, the decay length tends to -JD qp To. 
We also defined the dimcnsionlcss current Ii„j defined as 
the density of quasiparticles injected in the volume XA 
during the diffusion time A 2 / ' D qp . 

At low injection current, the above equations can be 
linearized. The excess in quasiparticle density is then di- 
rectly proportional to the injection current: z(0) f=s Ii n j. 
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FIG. 2: The dimensionless relative density of quasiparticles 
z(x) in logarithmic scale for two injection levels. The dimen- 
sionless current Ii„j is equal to respectively 0.1 (lower trace) 
and 1000 (upper trace). The horizontal axis is the x coordi- 
nate is in units of the diffusion length A. The a parameter is 
taken equal to 1. 



The density N qp decays exponentially with the distance 
from the interface, which gives a linear behavior on a log 
scale, see Fig. [2] lower trace. The decay characteristic 
decay length is the diffusion length A. 

At high current, i.e. when the injected quasiparti- 
clc density is larger than the thermal one, Eq. [7] gives 
z(0) « (3ali n j) 2 / 3 . This non-linearity originates from 
the recombination term N qp term in Eq. 2. The con- 
stant a is then small compared to 1 and the decay turns 
to a power law l/(x + a) 2 (see Fig. [2] upper trace). 
The crossover current between the two current regimes, 
i ss aDqpNqpoeA/ A, decreases with the thermal quasipar- 
ticle density. Therefore, the non-linearity persists even at 
low current when the temperature is low. 

Let us now discuss how the extra quasiparticles affect 
the operation of a N-l-S cooler. The presence of excess 
quasiparticles in the superconducting electrode has no 
direct signature on the current voltage characteristic of 
the tunnel junction. Essentially, the tunnel current is in- 
dependent of the quasiparticlc population in the super- 
conductor and can be written as an integral over energy 
involving only the normal metal distribution function f^: 
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Here n(E) is the BCS density of states ratio and is 
the normal state resistance of the cooler junction. Only 
under strong injection can a change in the I-V curve as- 
sociated with the gap suppression be observed. In con- 
trast, the heat current through the junction is an explicit 
function of the distribution function fs in the supercon- 
ducting side: 
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where both forward tunneling and back-tunneling pro- 
cesses contribute. 

In the usual analysis of the N-I-S junctions cooling 
power, the superconductor heating is ignored and there- 
fore the distribution function fs is taken to be equal to 
the Fermi function at the bath temperature To. In or- 
der to introduce the contribution of excess quasiparticles, 
one can include the actual non-equilibrium population in 
fs- As this population is peaked at the gap energy, it is 
sufficient to consider an additional heat current equal to 
A times the tunneling rate, so that: 
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SPm is the electronic thermal current due to the back- 
tunneling of the excess quasiparticles at x = 0. It is 
proportional to the excess quasiparticle density N qp (0) — 
N qp0 due to injection through the N-I-S junction. Besides 
the quasiparticlc back-tunneling, there is also a phonon 
contribution due to the reabsorption of 2A phonons in 
the normal electrode^ The excess 2A phonon density 
at position x = can be obtained as a function of N qp (0) 
from Eq. 2. It leads to a heat current ~ (N qp (Q) — N qp o) 2 . 
It is a second order process and will be ignored in the 
following. 

In a previous workjii we studied the electronic cooling 
of the central N-island in a S-I-N-I-S junction. The nor- 
mal island electronic temperature was extracted from the 
experimental current-voltage characteristics of the cooler 
junction. As expected, the experimental data show a 
minimum in electronic temperature at a bias just below 
the superconductor gap, when the heat current (Eq. 
is maximum. The cooling power was analyzed within a 
thermal model combining electron-phonon interactions in 
the normal metal and phonon transfer through the inter- 
face. At the optimum bias V « A/e, the electronic tem- 
perature was found to be systematically higher than pre- 
dicted. For example, at bath temperature of 292 mK, the 
lowest electronic temperature was about 97 mK whereas 
the prediction is 58 mK (see Fig. 4 of Ref. [TTh . We 
checked that the residual contributions due to imperfect 
filtering, leakage resistance or Andreev heating^ do not 
contribute significantly to the heating at the gap bias. 

Let us now examine how the heat losses due to the 
quasiparticle back-tunneling SP^ could explain the ob- 
served discrepancy. We estimate SPm from Eq. [TU] where 
the quasiparticle density has been substituted by its lin- 
ear approximation in Eq. [7] Within this approximation, 
the back-tunneling power is directly proportionnal to the 
tunnel current absolute value |i|: 
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The physical parameters of the device are a junction area 
A = 1.5 x 0.3 //m 2 , Cu and Al thickness dw = 50 nm and 
ds = 40 nm respectively, and the normal state resistance 
ryv = 1-25 kfi. The diffusion coefficient of the Al film 
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FIG. 3: Electronic temperature as a function of voltage bias 
in a S-I-N-I-S micro-cooler at a bath temperature of 292 mK. 
The data (circular dots) are fitted by the model described in 
Ref. [ll| (dotted line). The continuous line is obtained by 
adding the superconductor heating effect with the parameter 
/ = 14 pW/fiA. 

was measured at 4.2 K to be 30 cm 2 /s. The devices were 
designed so that every Al electrode is overlapped by a Cu 
strip, which plays the role of a normal metal trap. As a 
first analysis, we neglect its separation with the injection 
junction (0.3 fim) compared to the quasiparticle diffusion 
length (« 10 /an), so that we apply our simple analytical 
approach. The Al oxide barrier is identical in the trap 
and the cooler tunnel junctions, thus we can assume that 
their specific tunnel resistance is the same: Rtrap = rjyj4. 
Using N(E F ) = 3.44 x 10 10 eV~Vm~ 3 r^ and assum- 
ing that the Cu trap overlaps the superconductor over 
about half of its width, we find To = 0.4 fj,s. This is 
much smaller than the effective recombination time T e / / 
for aluminum^ It thus determines the effective diffusion 
length A w ^DgpTQ over which the quasiparticle energy 
is evacuated to the bath. If the free electron diffusion 
coefficient is used for D qp , we estimate the parameter / 



to be about 1 pW//LtA. 

We have added the contribution of SP^t = f\i\ to the 
thermal balance equation discussed in Ref. Illl using / as 
a fitting parameter and i equal to the experimental cur- 
rent. We found that for the bath temperatures 292 (see 
Fig. [3]) and 489 mK, the agreement with the experiment 
could be restored over the whole subgap bias voltage 
range by choosing / = 14 pW//iA. This is significantly 
larger than the calculated value. Still, we know that the 
latter is underevaluated for several reasons. First, we 
have overestimated D qp by taking the normal state value, 
which ignores the vanishing group velocity of quasipar- 
ticles near the gap energy. The reduction factor varies 
as feTjv/A at very low temperature. Since Pm is pro- 
portional to X/D qp , it must be multiplied by the inverse 
square root of this reduction factor, which amounts to 5 
at 0.1 A". Second, we have ignored the uncovered part of 
the superconducting electrode nearby the cooler junction. 
Eventually, alternative explanations of the efficiency mis- 
match at the gap edge remain fully possible. 

In summary, we developed a simple model for the dif- 
fusion of excess quasiparticles in a S-strip, including the 
effect of an external trap. Our assumption of a supercon- 
ducting line with an infinite length leads to simple and 
exact expressions. In more complicated designs used in 
experimental realizations, a numerical solution would be 
necessary. Our model can describe unexplained values of 
the minimum temperature observed in superconducting 
micro-coolers. It should also be useful to test accurately 
the influence of material and geometrical parameters and 
hence improve the performance of micro-coolers and su- 
perconducting tunnel junctions based low temperature 
detectors. 
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